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Abstract 



The inclusive y/i/r production cross-section and fraction of Jlijj mesons produced in B-hadron 
decays are measured in proton-proton collisions ai ^[s -1 TeV with the ATLAS detector at the 
LHC, as a function of the transverse momentum and rapidity of the JI41, using 2.3 pb ' of inte- 
. grated luminosity. The cross-section is measured from a minimum pj of 1 GeV to a maximum 

^ ' of 70 GeV and for rapidities within \y\ < 2.4 giving the widest reach of any measurement of 

I . production to date. The differential production cross-sections of prompt and non-prompt are 

separately determined and are compared to Colour Singlet NNLO*, Colour Evaporation Model, 
(— I . and FONLL predictions. 



1. Introduction 
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• The production of heavy quarkonium at hadron colliders provides particular challenges and 

\ opportunity for insight into the theory of Quantum Chromodynamics (QCD) as its mechanisms 

CO ■ of production operate at the boundary of the perturbative and non-perturbative regimes. Despite 

Tjlj- \ being among the most studied of the bound-quark systems, there is still no clear understanding 

■ of the mechanisms in the production of quarkonium states like the y/i/r that can consistently 

explain both the production cross-section and spin-alignment measurements in e^e^, heavy-ion 
and hadron-hadron collisions (see review articles [ JJ and references therein). 

Data obtained by the Large Hadron Collider (LHC) collaborations can help to test existing 
\r2 • theoretical models of both quarkonium production and /7-production in a new energy regime, 

\ at higher transverse momenta and in wider rapidity ranges than have previously been studied. 

d • Furthermore, quarkonium production in proton-proton collisions plays a key role as a reference 

point to understand heavy ion collisions and to understand the interplay between production and 
suppression mechanisms in such collisions 

This paper presents a measurement of the inclusive Jjij/ production cross-section and the 
production fraction fg of non-prompt Jlif/ (produced via the decay of a Z?-hadron) to inclusively- 
produced J 1 41 (hereafter referred to as the non-prompt fraction): 

_ cripp ^B + X^ J/if,X') 
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in the decay channel Jlijj IJ^/J^~ as a function of both y/i/r transverse momentum and rapidity in 
pp colHsions at the LHC at a centre-of-mass energy of 7 TeV and with an integrated luminosity 
of up to 2.3 pb The fraction has the advantage that acceptances and many efficiencies are 
the same for the numerator and denominator, and so systematic effects are reduced. The results 
of these analyses are compared to those made by the CMS Collaboration [3] with 314 nb ' of 
integrated luminosity and those from the CDF Collaboration |0] where appropriate. 

From these measurements, the prompt J/t// production cross-section {cr(pp — » J/t//X'), pro- 
duced directly from the proton-proton collisions or from decays of heavier charmonium states 
Uke the;ifc or t//(2S)), and the non-prompt {cr{pp B + X ^ J/ipX')) Jj^ production cross- 
section, are extracted. These results are compared to corresponding predictions made by the 
Colour Evaporation Model L5J, Fixed-Order Next-to-Leading Log (FONLL) [6J and Colour Sin- 
glet NNLO* calculations ITllsl]. Further details of the results of measurements presented here 
may be found in reference ||9|] . 



2. The ATLAS Detector and Data Processing 

In this section, the collection and processing of the data used in the paper are outlined. This 
involves a description of the most relevant subsystems of the ATLAS detector [10]: the trigger 
system, the muon system and the inner tracking detector Also specified are the triggers used and 
the offline data processing, in particular the selection of candidate muons. 

2.1. The ATLAS detector 

The ATLAS detector covers almost the full solid angle around the collision point with layers 
of tracking detectors, calorimeters and muon chambers. For the measurements presented in this 
paper, the trigger system, the inner detector tracking devices (ID) and the muon spectrometer 
(MS) are of particular importance. 

The ID covers the pseudorapidity range \ri\ < 2.5. It consists of a silicon pixel detector, a 
silicon strip detector (SCT) and a transition radiation tracker (TRT). These detectors are located 
at a radial distance from the beam axis between 50.5 mm and 1066 mm and are immersed in a 
2 T solenoidal magnetic field. The ID barrel consists of 3 pixel layers, 4 layers of double-sided 
silicon strip modules and 73 layers of TRT straws. The ID end-cap has 2x3 pixel layers, 2x9 
layers of silicon strips and 2 x 160 layers of TRT straws. 

The MS is located inside a toroidal magnetic field which provides 2.5 Tm of bending power in 
the barrel and 5 Tm in the end-caps. It consists of four detectors using different technologies and 
is divided into a barrel region (|?7| < 1.05) and two end-cap regions (1.05 < It/I < 2.7). Precise 
muon measurements are made using monitored drift tube chambers (MDT) in both the barrel 
and end-cap sections and using Cathode Strip Chambers (CSC) in the end-caps; fast triggers are 
obtained from resistive plate chambers (RFC) in the barrel and thin gap chambers (TGC) in the 
end-caps. The chambers are arranged in three layers, so high-p^ particles leave at least three 
measurement points with a lever arm of several metres. 

2.2. Trigger 

The ATLAS detector has a three-level trigger system: level 1 (LI), level 2 (L2) and the event 
filter (EF). For the measurements presented here, the trigger relies on the Minimum Bias Trigger 
Scintillators (MBTS) and the muon trigger chambers. 
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The MBTS are mounted in front of each liquid argon endcap calorimeter cryostat at z = +3.56 
m and are segmented into eight sectors in azimuth and two rings in pseudorapidity (2.09 < |?7| < 
2.82 and 2.82 < \ri\ < 3.84). The MBTS trigger is configured to require two hits above threshold 
from either side of the detector A dedicated muon trigger at the EF level is required to confirm 
the candidate events chosen for these measurements. This is initiated by the MBTS LI trigger 
and searches for the presence of at least one track in the entire MS. This trigger is referred to as 
the EF minimum bias trigger; it has an adjustable threshold on the reconstructed muon pr above 
which events are accepted and can be prescaled to accept a pre-determined fraction of events 
meeting the trigger condition. 

The LI muon trigger is based on RPCs for the barrel and TGCs for the end-caps [T^. It seeks 
hit coincidences within different RPC or TGC detector layers inside programmed geometrical 
windows which define the muon candidate pr, then selects candidates above six programmable 
thresholds and provides a rough estimate of their positions flT]. For the earlier data used in this 
analysis, the muon trigger corresponds to the lowest pr threshold trigger which requires a simple 
two-layer time coincidence within a region of 0. 1 x 0. 1 in ri-(p. No further geometrical constraint 
is applied. 

As the instantaneous luminosity of the collider increases, the trigger requirement switches 
from the EF minimum bias trigger to the LI muon trigger Later data periods make use of 
triggers seeded by this LI trigger but with additional pr cuts applied at the EF stage (these are 
referred to henceforth as the EF muon triggers). 

2.3. Muon identification and reconstruction 

Muon identification and reconstruction extends to I77I < 2.7, covering a pr range from 1 GeV 
up to more than 1 TeV. "Standalone MS tracks" are constructed entirely based on the signal hits 
collected in the MS. The track parameters are obtained from the MS track and are extrapolated 
to the interaction point, taking into account multiple scattering and energy loss in the traversed 
material. In this analysis, two categories of reconstructed muons are then defined: 

• Muons from combined reconstruction: the combined muon reconstruction relies on a 
statistical combination of both a standalone MS track and an ID track. Due to ID coverage, 
the combined reconstruction covers \ri\ < 2.5. 

• Muons from ID track tagging: a tagged muon is formed by MS track segments which 
are not formed into a complete MS track, but which are matched to ID tracks extrapolated 
to the MS. Such a reconstructed muon adopts the measured parameters of the associated 
ID track. In this paper, the muon tagging is limited to \ri\ < 2, in order to ensure high 
quality tracking and a reduction of fake muon candidates. 

The muon track helix parameters are taken from the ID measurement alone, since the MS does 
not add much to the precision in the lower momentum range relevant for the J/ifr measurements 
presented here. 

3. Data and Monte Carlo Samples 

Proton-proton collision data, at a centre-of-mass energy of 7 TeV, are included in this analysis 
if taken during stable beam periods and when the MS, ID and magnet systems were collecting 
data of a sufficiently high quality to be suitable for physics analysis. 
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Monte Carlo samples are used for determining acceptance corrections, as part of the trigger 
efficiency studies and in systematic cross-checks. They are generated using Pythia 6 |12] and 
tuned using the ATLAS MC09 tune 1 13] which uses the MRST LO* parton distribution functions 
1 14], The passage of the generated particles through the detector is simulated with Geant4 
OlSll and the data are fully reconstructed with the same software that is used to process the data 
from the detector. For the signal y/tA Monte Carlo (used to derive the kinematic acceptance 
corrections), the Pythia implementation of prompt y/i/r production sub-processes in the NRQCD 
Colour Octet Mechanism framework LI 6,1 is used. 

Prompt J 1x1/ production includes direct production from the hard interaction, as well as char- 
monium feed-down from excited states. These prompt production modes are distinct from non- 
prompt production that is characterised by the production of Jlxp via the decay of a B-hadron. 

All samples are generated with polar and azimuthal isotropy in the decay of the y/i/r (the de- 
fault in Pythia) and are reweighted at the particle level according to their respective angular de- 
pendencies in order to describe a number of different spin-alignment scenarios (see Section l4Tb . 
The J I if/ spin-alignment is not measured in this analysis, so the reweighted MC samples are used 
to provide an uncertainty band on the measurement of the production cross-section, determined 
by the maximum variation in acceptance across the full allowed range of Jj^ spin alignment. 

3.1. Event and candidate selection 

The analyses presented in this paper make use of the triggers described in Section l272l For 
the inclusive cross-section, in a given data taking period an event is retained or discarded based 
on the decision of a single specific trigger, without reference to any other triggers. For data from 
the initial running with lower instantaneous luminosity, the LI muon trigger is used. During later 
periods, with higher instantaneous luminosity, a more selective EF muon trigger with a 4 GeV 
Pt threshold is required, and eventually, this is increased to a 6 GeV pr threshold. The sample 
collected by these triggers and passing the data quality selections corresponds to an integrated 
luminosity of 2.2 pb"'. 

For the measurement of the B — > y/i/r non-prompt fraction (see Equation [TJ, two additional 
triggers are employed, and rather than using a single trigger to veto or accept events, several 
triggers are used simultaneously such that any one of them having fired results in the event 
being included. From the initial period, events triggering either the LI muon trigger or the 
EF minimum bias trigger are used (whereas only the LI muon trigger is used for the cross 
section). For intermediate instantaneous luminosities the LI muon trigger is used alone since 
the EF minimum bias trigger is highly prescaled at this stage. For the highest instantaneous 
luminosities, events are accepted which pass any of the EF muon triggers with pj thresholds of 
4, 6 or 10 GeV. During the runs with the highest instantaneous luminosities, the triggers with 4 
and 6 GeV are prescaled; however, the 10 GeV threshold trigger is not. The inclusion of this 
unprescaled trigger along with the addition of the EF minimum bias trigger for the B ^ J/iff 
non-prompt fraction measurement results in a slightly higher integrated luminosity of 2.3 pb '. 

To veto cosmic rays, events passing the trigger selection are required to have at least three 
tracks associated with the same reconstructed primary vertex. The three tracks must each have 
at least one hit in the pixel system and at least six hits in the SCT. 

Each remaining event is required to contain at least one pair of reconstructed muons. Only 
muons associated with ID tracks that have at least one hit in the pixels and six in the SCT 
are accepted. Di-muon pairs with opposite charges are considered to be J/iJ/ candidates if at 
least one combined muon is present in the pair At least one reconstructed muon candidate is 
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required to match a muon trigger (that is, at least one muon from the Jjijj candidate should 
have fired the trigger). For the early data, when the trigger is essentially based on the LI muon 
trigger, at least one of the offline muons is required to match the trigger muon candidate to within 
hR = ^fKif^^TK^ < 0.4 at the MS plane; for the later data taking, where the EF muon trigger is 
used, the offline and trigger muons are required to match within A/? < 0.005. 

The two ID tracks from each pair of muons passing these selections are fitted to a common 
vertex I17il. No constraints are applied in the fit and a very loose vertex quality requirement, 
which retains over 99% of the candidates, is used. 

For the B ^ J/tfr non-prompt fraction analysis, where lifetime information is an important 
element of the fit, additional requirements are made on the J/tp ^ A* A* candidates. The proba- 
bility of the fit to the J/i// vertex is required to be greater than 0.005. For this measurement J/i// 
candidates are rejected if the two muon candidate tracks were used to build different primary ver- 
tices in the offline reconstruction (so that there is an ambiguity as to which primary vertex to use 
in the lifetime calculation). This rejects fewer than 0.2% of the J /if/ candidates. This selection is 
not applied for the cross-section analysis. 



4. Inclusive J/i^ -* fi^fi Differential Production Cross-Section 

The measurement of the inclusive differential cross-section is determined as 

dMJI^) „ ,,,, + Nllt 

Bri JIip ^ fi^fi ) = — - — — (2) 
dpjdy -L ■ Apj-Ay 

where Nclf,. is the y/i/' yield in a given pr - y bin after continuum background subtraction and 
correction for detector efficiency, bin migration and acceptance effects, X. is the integrated lumi- 
nosity of the data sample and Apr and Ay are the pr and rapidity bin widths. The probability P 
that a J/if/ ^ fifi decay is reconstructed depends on the kinematics of the decay, as well as the 
muon reconstruction and trigger efficiencies. In order to recover the true number N^ifr of such 
decays produced in the collisions, a weight w is applied to each observed J/i// candidate, defined 
as the inverse of that probability and calculated as follows: 



-Jl-M- £l^^ ■ £%(p+, 7?^) ■ &-^{Pt, r]-) ■ £,rig (3) 



where ^ is the kinematic acceptance, Ai is a correction factor for bin migrations due to finite 
detector resolution, £trk is the ID tracking efficiency and fi^ is the single-muon offline reconstruc- 
tion efficiency. Here pj and 7/=^ are the transverse momenta and pseudorapidities of the positive 
and negative muons from the J /if/ decay. The trigger efficiency fitng for a given J/t// candidate is 
calculated from single-muon trigger efficiencies &f^g(Pj, J]"^) as follows: 

£tng = 1 - (l - StiM^ ■ (l - SmgC/'r^ V-)) ■ (4) 
The resultant weighted invariant mass peak is then fitted (see Section l4!4b to extract Ncifr- 
4.1. Acceptance 

The kinematic acceptance J[(pj,y) is the probability that the muons from a 7/ (A with trans- 
verse momentum pj and rapidity y fall into the fiducial volume of the detector This is calculated 
using generator-level Monte Carlo, applying cuts on the momenta and pseudorapidities of the 
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muons to emulate the detector geometry. Global cuts of |/?+|, |/?_| > 3 GeV for I77+I, [//-I < 2.5 are 
supplemented by finer pj thresholds in slices of 77 to ensure that regions of the detector where 
the values of offline and trigger efficiencies are so low as to be compatible with zero within the 
uncertainties (approximately 10%) are excluded from the analysis. 

The acceptance also depends on the spin-alignment of the y/t/r, which is not known for LHC 
conditions. The general angular distribution for the decay J/tfr ^ fifi in the J/i// decay frame is 
given by: 

oc 1 + cos^ 0* + Aj, sin^ 0* cos IS* + Agj, sin 26* cos S* (5) 

d cos 0*d(p* 

where 6* is the angle between the dkection of the positive muon momentum in the J /iff decay 
frame and the line of flight, while (p* is defined as the angle between the J/t// production and 
decay planes in the lab frame (see Figure [1] reference 111 811 and references therein). 

.Z* 

quarkonium 
rest frame 



production 
plane 




Figure 1: Definitions of the J/iJ/ spin-alignment angles, in the J/iJ/ decay frame. 8* is the angle between the direction of 
the positive muon in that frame and the direction of J/ip in the laboratory frame, which is directed along the z*-axis. (fi* 
is the angle between the J/ip production (x* - z*) plane and its decay plane formed by the direction of the J /if/ and the 
lepton f+ (from 1 18]). 

A large number of possible combinations of the coefficients Ag, A^, Ag^ have been studied, 
including some with Ag^ + 0. Five extreme cases have been identified that lead to the biggest 
variation of acceptance within the kinematics of the ATLAS detector and define an envelope in 
which the results may vary under all possible polarisation assumptions: 

1. Isotropic distribution, independent of Q* and 0*, with Ag - A^ - Ag^ = 0, labelled as 
"FLAT". This is used as the main (central) hypothesis. 

2. Full longitudinal alignment with Ag = -l,A^ = Ag^ = 0, labelled as "LONG". 

3. Transverse alignment with Ag = +1,A^ = Ag^ = 0, labelled as T+o. 

4. Transverse alignment with Ag = +l,A^ = +l,Ag^ - 0, labelled as T++. 

5. Transverse alignment with Ag = +l,A^ = -l,Ag^ - 0, labelled as T+_. 

Two-dimensional acceptance maps are produced in bins of pr and y of the for each of 
these five scenarios, and are illustrated in Figure |2] The maps are obtained by reweighting the 
flat distribution at the generator level using Equation |5] The central value for the cross-section 
measurement is obtained using the flat distribution, and the measurement is repeated using the 
other scenarios to provide an envelope of maximum variation, which is stated as a separate 
uncertainty. 
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Figure 2: Kinematic acceptance maps as a function of J/ip transverse momentum and rapidity for specific spin-alignment 
scenarios considered, wliicli are representative of tlie extrema of tlie variation of the measured cross-section due to spin- 
alignment configurations. Differences in acceptance beliaviour, particularly at low pr, occur between scenarios and can 
significantly influence the cross-section measurement in a given bin. 
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4.2. Bin migration corrections 

The measured efficiency and acceptance corrected Jjij/ pj distribution is parameterised in each 
rapidity slice by a smooth analytic function smeared with a Gaussian distribution, with resolution 
derived from the data. This function is integrated numerically over each analysis bin, both with 
and without smearing applied, and the ratio of the two integrals is assigned as the correction 
factor. The effects of this correction are minimal at low pj and at low rapidities (around 0.1%) 
but increase at higher pj and at higher rapidities (reflecting the decreasing momentum resolution) 
to a maximum of approximately 3%. 

4.3. Muon trigger and reconstruction efficiency 

The offline single muon reconstruction efficiencies are obtained from data using a tag and probe 
method 1.19.1 . where muons are paired with ID tracks ("probes") of opposite charge. The pairs 
are divided into two categories: those in which the probe is reconstructed as a muon ("matched") 
and those in which it is not ("unmatched"). Both sets of pairs are binned according to the pr 
and J] of the probe. In each of these bins, the muon reconstruction efficiency is obtained as 
the ratio of the number of J/tfr candidates in the peak of the matched distribution to the total 
number of candidates in the two mass distributions. The efficiency is extracted as a parameter of 
a simultaneous fit to both distributions. The dependence of the offline reconstruction efficiency 
on the muon charge is well described by MC within the acceptance. This procedure is repeated 
separately for combined and tagged muons. At higher pr (for muons with pr above 6 GeV), 
the efficiency determination is supported by additional tag and probe Z — > yu^/j" data L20i] for 
improved precision in the efficiency plateau region. 

A hybrid Monte Carlo and data-derived (tag and probe) scheme is used to provide trigger ef- 
ficiencies for the analysis with finer binning than would be possible with the available data statis- 
tics. This is necessary to avoid significant biases that would otherwise appear in the analysis with 
coarsely binned efficiencies across rapidly-changing efficiency regions. Due to significant charge 
dependence at low pr and high pseudorapidity, separate trigger efficiency maps are produced for 
positive and negative muons. Fully simulated samples of prompt pp — > J/i// (p^^pT ) X decays are 
used to populate the pr - y plane, using a fine binning. For each bin, the probability of a 
muon activating the trigger is determined. The derived efficiencies are then reweighted to match 
the data efficiencies in the reconstructed bins in cases where discrepancies exist between the data 
and Monte Carlo, and uncertainties from data are assigned. 

Muon reconstruction efficiencies have been determined relative to reconstructed ID tracks. 
Inner Detector tracks associated to muons and having the selection cuts used in this analysis have 
a reconstruction efficiency S^k of 99.5% + 0.5% per track (with no significant pseudorapidity or 
Pt dependence observed within the phase space probed with this analysis), which is applied as 
an additional correction to the J/tfr candidate yields. 

4.4. Fit of J lip candidate mass distributions 

The distribution of reconstructed J I if/ candidates over the candidate pj - y plane is shown in 
Figure |3] The majority of Jjili candidates are reconstructed in intermediate-pj-, high-y areas, as 
at lower pj values the acceptance of the detector is limited. 

The inclusive J I if/ production cross-section is determined in four slices of Jlijj rapidity: 
\y\ < 0.75,0.75 < \y\ < 1.5,1.5 < |y| < 2.0 and 2.0 < \y\ < 2.4. In Figure H the invariant 
mass distributions for all oppositely charged muon pairs passing the selection for the differential 
cross-section measurement are shown, before acceptance and efficiency corrections, for the four 
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Figure 3: Distribution of reconstructed Jjxp candidates (in the invariant mass interval 2.7 < niji^ < 3.5 GeV) as a 
function of J/if/ pj and rapidity. 



rapidity slices. Table [Upresents the results of the combined signal and background fits. In these 
fits the y/i/r and (/'(2S) peaks are represented by Gaussians, while the background is described by 
a quadratic polynomial. 



Table 1: Fitted mass, resolution and yields of Jjip candidates reconstructed in four Jjip rapidity bins. All uncertainties 
quoted are statistical only. The shift in mass away from the world average in the highest rapidity bin reflects the few-per- 
mille uncertainty in the tracking pj scale at the extreme ends of the detector. 





\y\ < 0.75 


7/ (A rapidity range 
0.75 < M < 1.5 1.5<M<2.0 


2.0 < \y\ < 24 


Signal yield 
Fitted mass (GeV) 
Fitted resolution (MeV) 


6710 ±90 
3.096 ± 0.001 
46+1 


10710+ 120 9630+ 130 
3.097 ± 0.001 3.097 + 0.001 
64+ 1 84 ± 1 


4130 ±90 
3.109 + 0.002 
111 + 2 
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Figure 4: Invariant mass distributions of reconstructed Jjifj —> n^n" candidates used in the cross-section analysis, 
corresponding to an integrated luminosity of 2.2 pb"' . The points are data, and the uncertainties indicated are statistical 
only. The solid lines are the result of the fit described in the text. 



The invariant mass distribution of Jjij/ — » ju^//" candidates in each pr - y bin is fitted using 
a binned minimum-;t'^ method. The J /if/ and t//(2S) signals are described by single Gaussians, 
while the background is treated as a straight line. 

For the differential cross-section measurement, the correction weight w defined in Equation[3] 
is applied to each candidate, and a new binned minimum-;i'^ fit is performed in each bin. The 
yields of J/i// determined from these fits, divided by the integrated luminosity, give the inclusive 
production cross-section for a given bin. Representative invariant mass distributions are shown 
in Figure |5] The probability distribution of the weighted fits across all bins is found to be 
consistent with the statistical expectation. 

The cross-sections obtained for each bin are listed in Table |2] the systematic uncertainties 
considered are displayed in Figure |6]and the cross-section results are presented in Figure Q The 
measurement in each pr - y analysis bin is positioned at the average pr for J/t// candidates in 
that bin. Various tests of the method described above are performed using simulated samples of 
known composition, and the number of J/t// in each analysis bin is successfully recovered within 
expectations in all cases. 

4.5. Systematic uncertainties 

Studies are performed to assess all relevant sources of systematic uncertainty on the mea- 
surement of the J/i// inclusive production cross-section. Sources of uncertainty are listed below, 
ordered according to the approximate size of their contribution (starting with the largest). 
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Figure 5: Acceptance- and efficiency-corrected invariant di-muon mass distributions scaled by integrated luminosity for 
selected bins in Jjtjj rapidity and transverse momentum. Low- and high- pj bins are shown here for the central and 
forward rapidity ranges, to represent the complete sample. Statistical uncertainties and systematic uncertainties due to 
efficiency and acceptance corrections are shown, combined in quadrature. 



1. Spin-alignment: Kinematic acceptance depends on the spin-alignment state of the 
and hence affects the corrected yield. Five spin-alignment scenarios are considered, which 
correspond to the extreme cases for the acceptance corrections within the kinematics ac- 
cessible in ATLAS. In each bin, the maximal deviations in either direction are assigned as 
the systematic uncertainty due to the unknown spin-alignment of the J I if/. These uncer- 
tainties are regarded as theoretical rather than experimental, and are quoted independently 
of the statistical and experimental systematic uncertainties. 

2. Muon reconstruction: The single muon efficiency maps are obtained from the data us- 
ing the tag and probe method, in bins of muon transverse momentum and pseudorapidity. 
Each efficiency has an uncertainty (predominantly statistical in nature, but with a system- 
atic component from the tag and probe method) associated with it. In order to obtain an 
estimate on the effects of uncertainties within these bins, the relative uncertainties (due to 
systematic and statistical components) on all Jjij/ candidates in a bin are averaged. Inner 
Detector tracks originating from muons and having the selection cuts used in this analysis 
have a reconstruction efficiency of 99.5% + 0.5% per track. The results are corrected for 
this efficiency, and a systematic uncertainty on the efficiency is assigned for each track, 
propagated linearly into the cross-section systematic. 

3. Trigger: The uncertainty on the trigger efficiency has components from the data-derived 
efficiency determination method (again largely statistical in nature) and from the reweight- 
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ing of MC maps to the data-driven (tag and probe) efficiency values. Tliese errors are 
treated similarly to the reconstruction efficiency uncertainties. 

4. Luminosity: The uncertainty on the integrated luminosity used for this measurement is 
determined to be 3.4% 112 111 , fully correlated between bins. 

5. Acceptance: 

• Monte Carlo statistics: The acceptance maps are obtained from dedicated Monte 
Carlo production, in bins of J 1 41 transverse momentum and rapidity. The acceptance 
in each bin has an uncertainty due to Monte Carlo statistics. The relative error on 
the acceptance correction for each y/f/r candidate contributing to a particular analysis 
bin is averaged in quadrature to evaluate the systematic effect of these errors on the 
cross-section measurement in that bin. 

• Kinematic dependence: The impact of any discrepancies in the underlying kinematic 
distribution modelling in the Monte Carlo used to build the maps, or differences in 
the pt dependence of the prompt and non-prompt components to the overall inclusive 
cross-section are studied. A correction to the acceptance maps is made based on 
the measured non-prompt to prompt fraction to ensure proper correction of the two 
populations, and an uncertainty is assigned based on the difference in yields from 
using the corrected and uncorrected maps. This uncertainty is significantly below 1 % 
in most analysis bins, reaching a maximum of 1.5% in some high pj, low rapidity 
bins. 

• Bin migration: The changes to the measured cross-section due to the migration of 
entries between the pj bins is determined by analytically smearing the efficiency and 
acceptance corrected pj spectrum with a Gaussian resolution function with width 
based on muon pj resolutions, taken from data. The correction needed to the central 
value due to bin migrations is as small as 0.1% at low pj and low rapidity and rises 
to ~ 3% at high pj and high rapidity. The variation of the bin migration correction 
within a given analysis bin (due to changing detector resolution and parameterisation 
of the Pt spectrum) is taken as a systematic. 

• Final-State Radiation: The acceptance maps correct the measured cross-section back 
to the y/iA kinematics, rather than the final-state muon kinematics, in order to allow 
proper comparison with theoretical predictions. Emission of QED final-state radia- 
tion is known to high accuracy, so the relative uncertainty on the modelling of this 
correction is determined to be less than 0. 1%. 

6. Fit: Invariant mass distributions for a large number of pseudo-experiments are constructed 
for each pr - y bin of the analysis, with the bin contents for each pseudo-experiment be- 
ing an independently Poisson-fluctuated value with mean equal to the measured data, and 
uncertainty in the bin determining the variance of the fluctuations. Within these pseudo- 
experiments, the candidate yields from the central fit procedure and yields from varied fit- 
ting models are determined, and the shift per pseudo-experiment calculated. The variation 
in fitting models include signal and background fitting functions and inclusion/exclusion of 
the tp{2S ) mass region. The means of the resultant shifts across all pseudo-experiments for 
each fit model are used to evaluate the systematic uncertainty. The largest mean variation 
in that bin is assigned as a systematic uncertainty due to the fit procedure. 

7. J/iJ/ vertex-finding efficiency: The loose vertex quality requirement retains over 99.9% 
of di-muon candidates used in the analysis, so any correction and systematics associated 
to the vertexing are neglected. 
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Figure 6: Summary of the contributions from various sources to the systematic uncertainty on the inclusive differential 
cross-section, in the Jlijj pj and rapidity bins of the analysis. The total systematic and statistical uncertainties ai'e also 
overlaid. The theoretical uncertainty due to the unknown spin alignment is not included on these plots. 

A summary of the various contributions to the systematic uncertainties on the measurement in 
each rapidity slice as a function of pj is shown in Figure |6] The uncertainty due to the 
luminosity (3.4%) is not shown, nor is the spin-alignment envelope which represents a full range 
of variation due to the unknown spin-alignment state. 

4.6. Inclusive J Itp cross-section results 

The results of the inclusive double-differential //(i_production cross-section measurement are 
given in Table|2] They are compared to CMS results fy\ in Figure|7]for cases where the rapidity 
ranges are close enough to permit comparison. The two sets of results show good agreement 
within experimental uncertainties and provide complementary measurements at low (CMS) and 
high (ATLAS) pt, together providing a measurement over a large kinematic range. 

The systematics are dominated by the data-driven muon reconstruction efficiency uncertain- 
ties, which are in turn dominated by their statistical uncertainties. There is an additional overall 
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uncertainty of +3.4% (fully correlated between bins) due to the luminosity measurement uncer- 
tainty. The measurement of the differential cross-section is limited by systematic uncertainties, 
with statistical uncertainties only contributing significantly near the low-pr thresholds where 
yields are limited by trigger efficiency, and in the highest transverse momentum bin. 

The total cross-section for inclusive J/i// — > /u^fi^ production, multiphed by the branching 
fraction into muons and under the FLAT production scenario for the central value, has been 
measured for J /iff produced within |y| < 2.4 and pj > 1 GeV to be: 

Br(7/iA H^H-)(T(pp J/iffX; \yj/^\ < 2A, > 7 GeV) 

= 81 ± 1 (stat.) ± 10(syst.) ±11 (spin) ± 3 (lumi.) nb 

and for J/t// within 1.5 < lyl < 2 and pr > I GeV to be: 

Br(J/il^ ^ ix'n-)(T{pp ^ JlijfX- 1.5 < \yji^\ < 2,pf > 1 GeV) 

= 510 + 70 (stat.) +1% (syst.) +lf^ (spin) ± 20 (lumi.) nb. 



5. Measurement of the Non-Prompt J /if/ Fraction 

Experimentally, it is possible to distinguish J/ ^ from prompt production and decays of heav- 
ier charmonium states from the Jlijj produced in B-hadron decays (non-prompt production). The 
prompt decays occur very close to the primary vertex of the parent proton-proton collision, while 
many of the //i^ mesons produced in Z?-hadron decays will have a measurably displaced decay 
point due to the long lifetime of their Z?-hadron parent. 

From the measured distances between the primary vertices and corresponding Jl^ decay 
vertices the fraction fg of Jl^ that originate from non-prompt sources, as defined in Equation[T] 
can be inferred. An unbinned maximum likelihood fit is used to extract this fraction from the 
data. 



5.7. Pseudo-proper time 

The signed projectioi 
constructed according to the following formula 



The signed projection of the Jjip flight distance, L, onto its transverse momentum, p'j'^ , is 



L,,,^L-/4^^/pf, (6) 

where L is the vector from the primary vertex to the Jltp decay vertex and p'j'^ is the transverse 
momentum vector of the Here Lw measures the displacement of the Jl^ vertex in the 

transverse plane. 

The probability for the decay of a B-hadron as a function of proper decay time t follows an 
exponential distribution 

p{t) = — exp(-f/rB), (7) 

where tb is the lifetime of the B-hadron. For each decay the proper decay time can be calculated 
as 

L 

t = (8) 
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Figure 7: Inclusive J/if/ production cross-section as a function of J/i// transverse momentum in the four rapidity bins. Overlaid is a band representing the variation of the result 
under various spin-alignment scenarios (see text) representing a theoretical uncertainty. The equivalent results from CMS are overlaid. The luminosity uncertainty (3.4%) is 
not shown. 



where L is the distance between the B-hadron production and decay point and jiy is the Lorentz 
factor. Taking the projection of the decay length and momentum on the transverse plane for 
B-hadrons, one obtains 



Pt 



(9) 



In this case, L^y is measured between the position of the reconstructed secondary vertex and the 
primary vertex in the event. The primary vertex is refitted with the two muon tracks excluded, 
to avoid a bias. The uncertainty on L^y is calculated from the covariance matrices of the primary 
and the secondary vertices. The majority of the events contain only a single primary vertex. In 
the few that contain multiple vertices, the J /if/ is assigned to a primary vertex based on the use of 
the tracks by the ATLAS reconstruction software; if both J/t// tracks are included in the recon- 
struction of the same primary vertex, this is the one which is assigned. In a small number of cases 
(fewer than 0.2%) the two tracks making the J/i// candidate are included in the reconstruction of 
different primary vertices. These candidates are discarded. 

Since the B-hadron is not reconstructed completely, one does not know its transverse mo- 
mentum. Instead the J/il/ momentum is used to construct a variable called the "pseudo-proper 
time" 

r=^. (10) 
Pt 

Here, the world average value of ni'p^Q is used to reduce the correlation between the fits that will 
be performed on the mass and the Ufetime. Studies show that the results are insensitive to this 
choice. 

At large pj'^, where most of the B-hadron transverse momentum is carried by the J/ifr, the 
distribution of t is approximately exponential, with the B-hadron lifetime as a parameter. At 
small Pj'^, the range of opening angles between the J/i// and B-hadron momentum leads to a 
smearing of the underlying exponential distribution. 

5.2. Fitting procedure 

The sample is divided into bins of pt and rapidity y of the Jjij/ candidates. In each bin, 
a maximum likelihood fit is performed in order to determine the fraction of the non-prompt to 
inclusive Jlijj production cross-sections in that particular bin. The mass and pseudo-proper time 
are simultaneously fitted in the entire mass region from 2.5 to 3.5 GeV, using the likelihood 
function: 

N 

L = W [/sig^'sig(T, Sr)Fsigim^i„ (5m) + (1 - /sig)Pbkg(T, 5r)f'bkg(w^//)] (1 1) 

where is the total number of events in the 2.5 - 3.5 GeV mass region and /sig is the fraction 
of signal J/il/ candidates in this region determined from the fit. Pgig and Pbkg are pseudo-proper 
time probability density distributions (PDFs) for the J/i// signal and background candidates re- 
spectively, and are described fully below. The F^ig, Fbkg functions are the mass distribution 
models for signal and background. In sunmiary, the input variables to the maximum UkeUhood 
fit to determine the production ratio are the pseudo-proper time t, its uncertainty 6t, the di-muon 
mass m^^ and its uncertainty 6m. 
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5.2.7. Invariant mass and pseudo-proper time probability density functions 
For the signal, the mass is modelled with a Gaussian distribution: 

Fsig(m^^i,6„,)= —= e -(^''."'^ (12) 

yln S6,„ 

whose mean value nij/^ is the J/i// mass, determined in the fit, and whose width is the product 
SSm, where 6,„ is the measured mass error calculated for each muon pair from the covariance 
matrix of the vertex reconstruction and 5 is a global scale factor to account for a difference 
between 6,„ and the mass resolution from the fit. For the background, the mass distribution is 
assumed to follow a second-order polynomial function. 

The pseudo-proper time PDF for J/ifr signal candidates, f sig, consists of two terms. One 
term describes the J/iJ/ from B-hadron decays (Pb), and the other describes the J/i// from prompt 
decays (Pp): 

P.ig(T, Sr) = fsPBiT, Sr) + (l- fB)Pp{r, <5.), (13) 

where fs is the fraction of Jjijj from B-hadron decays as defined in Equation[T] 

The pseudo-proper time distribution of the 7/ ^ particles from B-hadron decays Pb{t, 5t) is an 
exponential function E{t) - exp(-T/Tett) with a pseudo-proper time slope Teff, convolved with 
the pseudo-proper time resolution function R{t' - t, 5^): 

Pb(t, 6r) = RiT - T, 6r) ® E(t'). (14) 

Promptly produced J/ifr particles decay at the primary vertex, and their pseudo-proper time dis- 
tribution is thus given by the resolution function: 

Pp(T, 6r) = R(t' - T, 6r) ® 6(t') = RiT, 6r). (15) 

The resolution function /? is a Gaussian distribution centred at r = with a width S ,6t, where S , 
is a scale factor (a parameter of the fit) and dj is the per-candidate uncertainty on t, the measured 
pseudo-proper lifetime determined from the covariant error matrix of the tracks. 

The pseudo-proper time PDF for background candidates Pbk? consists of the sum of a long 
lived component modeled with an exponential function and a prompt component modeled by 
a delta function and two symmetric exponential tails. Each component is convolved with the 
Gaussian resolution function: 

Phko{T,dr) = J(l -bi- b2)6{T') + bi expj;^j + /?2exp J^jj ®/?bkg(T' - T,6r), (16) 

where Rbkg(T) is a Gaussian distribution centered at t = with a width S^kg^T, where Sbkg is a 
scale factor (a parameter of the fit) and Sj is the per-candidate uncertainty on the measured t. 
Parameters Teffi and are pseudo-proper time slopes of the two components of background, 
and bi and b2 are the corresponding fractions of the background. All four parameters (Teffi, Teif2, 
bi and bi) are determined from the fit. 

5.2.2. Summary of free parameters 

The full list of the parameters of the fit are as follows: 

• /sig the fraction of signal J/t// candidates in the 2.5-3.5 GeV mass region of the fit; m j/^ 
the mean value of the J/t// mass; the scale factor S to account for a difference between 5„, 
and the mass resolution from the fit; 
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• fs the fraction of from B-hadron decays; a pseudo-proper time slope Teff describing the 
B-hadron decays; S , a scale factor to account for a difference between 6t and the Z?-hadron 
pseudo-proper time resolution from the fit; 

• the slope parameters Teifi, Teif2 and 5bkg describing the time evolution of the Jlijj back- 
ground, in analogy to the parameters of B-hadron decays, defined above; bi and b2, frac- 
tions of the two background components. 

5.3. Results of the likelihood fits 

The results of the likelihood fit to the pseudo-proper time distributions in a representative 
Pj'^ bin are shown in Figure|8] The figure shows the result of the unbinned maximum likelihood 
fits for the signal and background components projected onto the lifetime and invariant mass 
distributions. From the results of the fit, it is possible to derive the non-prompt to inclusive 
production fraction as a function of Pj'^. The probabilities and Kolmogorov-Smimov test 
results for the fits across all analysis bins are found to be consistent with statistical expectations, 
with the lowest fit probability out of over 70 fits being 1%. 
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Figure 8: Pseudo-proper time distributions (left) of Jjip — > ^i^fi" candidates in the signal region, for a selected pj bin 
9.5 < Pj < 10.0 GeV in the most central and most forward rapidity regions. The points with error bars are data. The 
solid line is the result of the maximum likelihood unbinned fit to all di-muon pairs in the 2.5 - 3.5 GeV mass region 
projected on the naiTow mass window 2.9 - 3.3 GeV. The invariant mass distributions which are simultaneously fitted 
with the pseudo-proper time are shown on the right for the same bins. 
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5.4. Systematic uncertainties 

Several studies performed to assess all relevant sources of systematic uncertainties on the 
measured fraction of non-prompt to inclusive J/tJ/ decays are outlined below, in order of impor- 
tance. 

1 . Spin-alignment of prompt J/t//: In general, spin-alignment may be different for prompt 
and non-prompt J/i//, which may result in different acceptances in the two cases. The 
central value assumes they are the same (isotropic distribution in both angles, as for the 
inclusive cross-section central result), but four additional scenarios for the prompt com- 
ponent are also considered, as discussed in Section 14.11 The largest variations within the 
four models from FLAT is calculated for each bin in turn and assigned as an uncertainty 
envelope on prompt production. 

2. Spin-alignment of non-prompt J/t//: The possible variation of spin-alignment in B — » 
J/t//X decays is expected to be much smaller than for prompt J/i// due to the averaging 
effect caused by the admixture of various exclusive B — > J/i/rX decays. We assign an 
additional uncertainty on the non-prompt fraction (and non-prompt cross-section) for the 
difference in final result when using either an isotropic spin-alignment assumption for non- 
prompt decays or maps reweighted to the CDF result ll22t for B ^ J/i/i spin-alignment. 
This contributes up to an additional 0.4% uncertainty on the overall (prompt and non- 
prompt) systematic due to spin-alignment on the fraction. 

3. Fit: A number of changes are applied to the fitting procedure, and the fit is repeated in 
order to gauge the sensitivity of the fraction to the details of the fits: 

• The central value for the fraction assumes a background model for the proper time 
distribution of the background that includes one exponential function with a negative 
slope and a symmetric double exponential term with the same absolute value, Teif2, 
for the negative and positive slopes. To test the robustness of the result, this model 
is changed in two ways. First, the symmetric term is no longer required to be sym- 
metric, so different values of the negative and positive slopes are allowed. Second, 
the sum of two asymmetric double exponentials is used, having the same negative 
decay constant but differing positive decay constants. The maximum deviation from 
the central value is taken as a systematic uncertainty. 

• The per-candidate Gaussian convolution function is changed to a per-candidate dou- 
ble Gaussian convolution, allowing different scale factors (to account for differences 
between the resolution returned by the tracking algorithm and measured resolution) 
for each Gaussian to be determined from the fit. Differences from the main fit are 
assigned as a systematic uncertainty. 

• The main result uses a second-order polynomial in the mass fit to describe the back- 
ground. To test the sensitivity to this choice, the fits are repeated using instead poly- 
nomials of degree one and three. Differences from the main fit are assigned as a 
systematic. 

• The central result takes J/t// candidates in a mass range from 2.5 to 3.5 GeV, to avoid 
the mass region of the i/'(2S). In order to test the stability of the result and to increase 
the statistics in the side bands, the analysis is repeated with a mass range from 2 to 
4 GeV, but excluding the region from 3.5 to 3.8 GeV. The result is stable compared 
to the statistical uncertainties, and so no systematic uncertainty is assigned for this 
source. 
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• The analysis relies on a simultaneous fit to the proper time and mass distributions. 
The likelihood used assumes no correlation between the two quantities. To test the re- 
liability of this assumption, the mean measured invariant mass is plotted as a function 
of the proper time. The resulting distribution is flat, except in the negative lifetime 
region and at very long proper lifetimes, where residual background dominates the 
sample and invalidates the test. Accordingly, no expUcit systematic for this correla- 
tion is assigned. 

4. Kinematic dependence: Differences in the acceptance of prompt and non-prompt Jjij/ 
due to their different momentum spectra, averaged across an analysis bin, can bias the 
fraction measurement. A correction factor is calculated based on the acceptance maps 
with and without momentum reweighting to account for the differences between prompt 
and non-prompt J I if/ and this correction assigned as a systematic uncertainty. 

5. Reconstruction efficiencies: The central result for the fraction assumes that the recon- 
struction efficiencies are the same for non-prompt and prompt y/i/r mesons and hence can- 
cel in extracting the fraction. This assumption is tested on Monte Carlo samples described 
in Section [3] and no statistically significant shift is observed. Thus, no systematic uncer- 
tainty is assigned. 

6. PUe-up/multiple interactions: Some collisions result in the reconstruction of multiple 
primary vertices. The primary vertex chosen determines the transverse decay displacement 
Lxy used in the proper time determination. The central value is obtained by taking the 
primary vertex that is formed using both of the Jlijj candidate muons and rejecting cases 
where those candidates are associated with different primary vertices. To assess the effect 
of this procedure, two alternate methods where used. The first chooses the primary vertex 
with the highest summed squared transverse momenta of the tracks that form it. The 
second takes the same vertex, but rejects cases where either of the muon candidates are not 
used in determining that primary vertex. As no significant variation is seen in the results 
from the two methods, no additional uncertainties are assigned due to this source. 

The stability of the method used is checked using simplified Monte Carlo trial experiment 
samples to perform various tests of the closure of the analysis. The simultaneous mass and 
pseudo-proper time fit model is used to generate 100 simplified Monte Carlo experiments for 
each pj and y bin. The number of events generated is approximately the same as the number of 
data events for the corresponding bin. For each event the invariant mass and pseudo-proper time 
values are generated randomly from the total PDF, while the per-candidate error on invariant mass 
and pseudo-proper time are sampled from the corresponding experimental data distributions. 

For each experiment, a fit of the total PDF on the simple Monte Carlo sample is performed. 
The pull. A, defined as 

^_ (/generated -/extracted) 
0"(/extracted) 

is computed for each Monte Carlo experiment. Here /generated is the non-prompt fraction for the 
signal component according to which the Monte Carlo samples are generated (i.e. the result of 
the fit of the global model to the experimental data), while /extracted and cr(/extracted) are the value 
and uncertainty obtained from the fit. The Gaussian mean and sigma are statistically compatible 
with zero and unity, respectively, in all bins, indicating that no bias or improper uncertainty 
estimate is introduced by the fit. 
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5.5. Fraction of non-prompt Jjxli as a Junction of J/i// transverse momentum and rapidity 

Figure |9] and Tables[3]to|6]show the results of the differential non-prompt fraction measure- 
ment as a function of average p^'^, in each of the four rapidity bins. The uncertainty envelopes 
due to the unknown spin-alignment are overlaid as solid bands. 

The measurements are compared with those of CMS [3] and CDF [4] and build upon those 
results with finer rapidity binning, a much extended rapidity coverage relative to CDF and sig- 
nificantly increased pr reach relative to both experiments. Strong pr dependence of the fraction 
is observed: ~ 90% of J /if/ are produced promptly at low pr, but the fraction of non-prompt J/t// 
rapidly increases at mid-pr from ~ 15% at 7 GeV to ~ 70% at the highest accessible pr values. 
No significant rapidity dependence is seen. The ATLAS results exhibit good agreement with 
CMS results where they overlap, and also with the CDF measurements, indicating that there is 
no strong dependence of the fraction on collision energies. 

6. The Prompt and Non-Prompt Differential Production Cross-Sections 

The prompt and non-prompt J/tfr production cross-sections can be derived from the inclusive 
production cross-section and the non-prompt fraction. Where necessary, pr bins in the inclusive 
cross-section are merged to align bins in the prompt/non-prompt cross-section result with those 
in the non-prompt fraction measurement. The relative systematic uncertainties in each of the 
fraction and inclusive cross-section measurement bins (merged where appropriate) are taken to 
be uncorrected, while the statistical uncertainties are combined taking correlations into account. 
The spin alignment uncertainties are quoted independently of the experimental uncertainties. 

6.1. Non-prompt differential production cross-sections 

We assume the spin-alignment of a J/ifr meson from a B — > J/tfrX decay has no net polar or 
azimuthal anisotropy for the central result, as the possible variation of spin-alignment in Z? — > 
J/t/^X decays is expected to be much smaller than for prompt due to the averaging effect 
caused by the admixture of various exclusive B — > J/tfrX decays. We assign a spin-alignment 
uncertainty on the non-prompt cross-section for the difference in the final result when using either 
an isotropic spin-alignment assumption for non-prompt decays or maps reweighted to the CDF 
result 1 22] for B — > J/t// spin-alignment. 

The total integrated cross-section for non-prompt J/i/r, multiplied by the branching fraction 
into muons and under the "FLAT" production scenario, has been measured for J/t// mesons pro- 
duced within \y\ < 2.4 and pr > 7 GeV to be: 

Br(J/t// fi^p-)o-{pp B + X ^ J/t/jX; \yji^\ < 2.4, p^'^ > 7 GeV) 

= 23.0 + 0.6 (stat.) ± 2.8(syst.) ± 0.2 (spin) ± 0.8 (lumi.) nb 

and for J/t/i mesons produced with L5 < |y| < 2 and > 1 GeV to be: 

Br(J/^ p^ti-)o-ipp ^ B + X^ J/t//X; 1.5 < \yji^\ < 2,/?^^^ > 1 GeV) 

= 61+24 (stat.) + 19 (syst.) + 1 (spin) ± 2 (lumi.) nb. 



21 



g 
o 

"O 

o 



0.9^ 
0.8^ 
0.7^ 
0.6^ 
0.5^ 
0.4^ 
0.3^ 
0.2^ 
0.1^ 
0^ 



~n ' ' ' 

, ATLAS \/s=7 TeV, |yj,^l<0.75 
□ CMS\/s=7TeV, |y |<1.2 
o CDF \fs=1 .96 TeV, |yjj<0.6 
^Spin-alignment enveiope 



ATLAS 

f L dt ~ 2.3 pb ' 



10 



pf [GeV] 



c 1 

g 

0.9 

nJ 

^ 0.8 
o 

0.7 

"§ 0.6 

1 0.5 
^ 0.4 

CL 

E 0.3 

I 0.2 

i 0.1 




, ATLAS \/s=7 TeV, 0.75<ly |<1.5 

J/\|/ 

□ CIVIS \/s=7 TeV, |y |<1.2 

J/\|/ 

o CIVIS \ls=7 TeV, 1 ■2<|yj^^|<1 .6 
Spin-aiignment envelope 



ATLAS 

jLdt~2.3 pb ' 



10 



pf [GeV] 



c 
o 

o 
S? 

c 
g 

"O 

o 



Q. 

E 
o 



Ip- 
0.9^ 
0.8^ 
0.7^ 
0.6^ 
0.5^ 
0.4^ 
0.3^ 

0-2 h 
0.1^ 

0^ 



T" 



T" 



, ATLAS \/s=7 TeV, 1.5<|y |<2.0 
□ CMS\/s=7TeV, 1.6<|yj^^J<'2.4 
Spin-alignment enveiope 



ATLAS 

f L dt ~ 2.3 pb ' 



^1 



10 



pf [GeV] 



g 

o 
cd 

g 

o 

"a 
o 



1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




, ATLAS n/s=7 TeV, 2.0<|y 1<2.4 
□ CIVIS \/s=7 TeV, 1 .6<|yj^^|<2.4 
Spin-aiignment envelope 



ATLAS 

jLdt~2.3 pb ' 



^ 




10 



pf [GeV] 



Figure 9: J/i// non-prompt to inclusive fractions as a function of J/i// transverse momentum. Overlaid is a band representing the variation of the result under various spin- 
alignment scenarios (see text) representing a theoretical uncertainty on the prompt and non-prompt J/t// components. The equivalent results from CMS yj and CDF ||4j| are 
included. 



6.7.7. Comparisons with theoretical predictions 

ATLAS non-prompt J/tfr production cross-section measurements are compared to Fixed Or- 
der Next-to-Leading Logarithm (FONLL) calculations [6] in Tables |7] and [TO]and in Figure [TOl 
FONLL vL3.2 is used for these predictions, using the CTEQ6.6|23] parton density function 
set. FONLL predictions use a B — » J/iJ/X branching fraction of Br(B — » J/tfr) = 0.0116. Un- 
certainty bands associated with the predictions come from the input b-quark mass, varied within 
4.75 ±0.25 GeV, renormalisation and factorisation (fif) scales (independently) varied within 
0.5 < pR,f/m < 2 (with the additional constraint that 0.5 < UrI^f < 2) and parton density func- 
tion uncertainties. Good agreement is seen between the experimental data and the theoretical 
prediction across the full range of rapidity and transverse momentum considered. 

6.2. Prompt differential production cross-sections 

The prompt production cross-section is of direct interest for the study of quarkonium pro- 
duction in QCD. The spin-alignment state and pr dependence of the spin-alignment of promptly 
produced J/il/ particles are thought to be non-trivial, so the spin-alignment uncertainty envelope 
on the inclusive cross-section measurement is propagated into the prompt cross-section measure- 
ment. The prompt production cross-sections are presented in Tables fTTItofT?! 

The total cross-section for prompt J/tfr (times branching fraction into muons) under the flat 
production scenario has been measured for J/i// produced within |y| < 2.4 and pr > 7 GeV to be: 

Br(J/ifr n^ijr)a-{pp prompt J/i//X; \y\ < 2.4, pr > 7 GeV) 

= 59+1 (stat.) + 8(syst.) + I (spin) ± 2 (lumi.) nb 

and for J/i/r within 1 .5 < lyl < 2 and pr > I GeV to be: 

Br(J/i/r fi^jj-)o-(pp prompt J/ifrX; 1.5 < \y\ <2,pT>l GeV) 

= 450 ± 70 (stat.) , (syst.) ± (spin) + 20 (lumi.) nb. 

6.2.7. Comparisons with theoretical predictions 

In Figure[TT]the prom pt p roduction data are compared to the predictions of the Colour Evap- 
oration Model (CEM) lf5i l24ll for prompt J /if/ production (with no uncertainties defined) and a 
calculation of the direct production cross-section in the Colour Singlet Model (CSM) |25, 2^ 
at next-to-leading order (NLO) and a partial next-to-next-leading order calculation (NNLO*). 

The Colour Evaporation Model predictions are produced using the CTEQ6M parton density 
functions, a charm quark mass of 1.2 GeV and the renormalisation and factorisation scales set to 

= 2 ^p^ + TOg + kj (where pr is the transverse momentum of the J/i// and mg is the quark 

mass and kj is a phenomenological fit parameter set to 1 .5 GeV^). The CEM predictions include 
contributions from;t'c and i^(2S) feed-down and can be directly compared with the prompt J/tf/ 
data. The normalisation of the CEM prediction is generally lower than in data and strongly 
diverges in shape from measured data, showing significant disagreement in the extended pr 
range probed by the measurement described in this paper 

The Colour Sing let NLO and NNLO* prediction^] for direct J/tfr production use a charm 
quark mass of 1.5 GeV, the CTEQ6M parton density function set, and factorisation and renor- 



'Xhe NNLO* calculation is not a full next-to-next-to-leading order prediction, as it does not include all loop correc- 
tions to p/) — > Q + jjj (where j is a light parton) up to order a^. This limits the applicability of the calculation to values 
above a particular J/i/r pr threshold (due to soft and coUinear divergences). 
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malisation scales set to /^o = ^Pj + (varied up and down by a factor of two to determine 
scale uncertainties). As the calculation is for direct production, corrections must be applied for 
Xc and ip{2S) feed-down to bring the calculations in direct comparison with data. To correct for 
feed-down, a flat 10% correction is applied to account for the contribution of i/'(2S) — > J/tfrnn 
and a 40% correction is added to account for radiative Xc decays. This yields a total correction 
of 50%. The correction factor is not well-determined from theory or experiment so is assigned a 
100% uncertainty. This uncertainty is not included in the CSM theoretical uncertainty. 

The NLO and NNLO* predictions are overlaid with the ATLAS measurements in FigurefTTI 
for each rapidity region. The dashed lines represent the central NLO and NNLO* predictions 
while the shaded areas show the range of the prediction due to factorisation and renormalisation 
scale variation (although the upper band of this uncertainty may not encapsulate the full range of 
infrared uncertainties [7|). 

The Colour Singlet Model predictions at NNLO* show significant improvement in describing 
the pt dependence and normalisation of prompt production over NLO, and vast improve- 
ment over earlier LO predictions that are compared to Tevatron data, although it is clear that 
these predictions still fall short of fully describing the production mechanisms of prompt J /if/, 
particularly at the highest transverse momenta explored in this analysis. The overall scale of the 
central prediction is somewhat low, but these discrepancies are similar in nature to those seen 
between NNLO* calculations and i/'(2S) production as measured by CDF ||26[ I27I1 at lower pj 
and centre-of-mass energy and may be attributed to higher order corrections beyond NNLO* that 
are still expected to be relatively significant for hidden charm production. 

7. Summary 

Results are reported on the measurement of the inclusive cross-section of J/t// p^p^ pro- 
duction in proton-proton collisions at a collision energy of 7 TeV using the ATLAS detector with 
up to 2.3 pb"' of integrated luminosity. The inclusive cross-section is measured in bins of rapid- 
ity y and transverse momentum pr of J /if/, covering the range \y\ < 2.4 and 1 < p/- < 70 GeV. 
The fraction of non-prompt J/i// mesons is also measured as a function of J/i// transverse mo- 
mentum and rapidity and using the above two measurements, double-differential cross-sections 
are extracted separately for promptly-produced J/i// mesons and those coming from Z?-hadron 
decays. 

It is found that the measurements made by ATLAS and CMS are in good agreement with 
each other in the overlapping range of moderate pr values and complement each other at high 
(ATLAS) and low (CMS) values of transverse momenta. The non-prompt production fraction 
results are also compared to those from CDF at lower energy and reasonable agreement is found, 
suggesting there is no strong dependence of the fraction on the collision energy. 

The results are also compared to various theoretical calculations of prompt as well as non- 
prompt J/i// production. In general, the theoretical curves describe the non-prompt data well, 
but significant deviations are observed in the prompt production spectra both in shape and nor- 
malisation, particularly at high transverse momenta. These measurements can thus provide input 
towards an improved understanding and theoretical description of J /if/ hadronic production. 
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Table 2: Inclusive Jj^ji production cross-sections as a function of Jltfi pr in four rapidity (|y|) bins. The first uncertainty 
is statistical, the second is systematic and the third encapsulates any possible variation due to spin-alignment from the 
unpolarised (ie = = Ag^ = 0) central value. 



(Pt) 
(GeV) 



Value 



2 < |>.| < 2.4 
± (stat.) ± (syst.) 



dpjdy 

± (spin) 



Br( Jlip P+H-) [pb/GeV] 



(pt) 




1.5 


(GeV) 


Value 


± (stat.) 


2.8 


143000 


±23000 


4.5 


39400 


±5500 


5.2 


17600 


±3300 


5.7 


14300 


±1200 


6.3 


12760 


±920 


6.8 


8910 


±610 


7.2 


6350 


±430 


7.7 


5040 


±350 


8.3 


3790 


±210 


8.7 


3110 


±160 


9.2 


2260 


±110 


9.7 


1674 


±85 


10.5 


1297 


±46 


11.5 


754 


±31 


12.9 


404 


±15 


14.9 


193 


±10 


16.9 


103.0 


±6.9 


19.6 


48.9 


±3.2 


25.0 


10.6 


±1.1 


34.1 


2.22 


±0.40 



■■\y\<2 

± (syst.) 
rism — 



± (spin) 
274000 



5.3 
5.8 

6.3 
6.8 
7.3 
7.7 
8.3 
8.7 
9.2 
9.7 
10.5 
11.5 
12.9 
14.9 
16.9 
19.7 
24.9 
33.6 



15900 
13500 
8800 

6290 
3990 
4070 
2650 
1930 
1450 
1208 
829 
598 
320 
164 
77.8 
29.9 
6.2 
1.12 



±4300 
±3600 
±1100 

±700 

±500 

±450 

±290 

±160 

±130 

±94 

±51 

±43 

±19 

±12 

±8.2 

±3.3 

±1.1 

±0.43 



.2800 



-550 
,570 



-390 
,260 



-138 

-92 
+69 
-73 
+38 

-8.0 
±3-1 

.3:6 
% 

0.7,8 



,28800 

-2200 
,5140 

-m 

-690 
,2920 



-166 

,286 

-87 

^.174 

-71 

-40 

1-33 
-21 
^14.1 



8:06 

10 



M 

860 
,590 



M 

183 
,146 
-139 

+90 
-83 
+« 
-43 

+21 
-19 
,13.0 
-9.4 
^4.1 
^4.2 
^1.0 

M 
0.7,l 



M 

1430 
.2260 

1% 



M 

296 
,316 
-241 
,168 
-147 

-75 
1-28 
-32 
,12.0 



±1:2 
.0.13 



(pt) 
(GeV) 



Value 



0.75 < 
± (stat.) 



\y\ < 1.5 
± (syst.) 



,4100 
-3800 
,2700 



-1040 

,860 

-810 

,420 

-440 

,390 

-390 

.330 

-320 

^.i50 

;m 

;m 
-120 

,60 



.0.34 

'Ml 

-0.043 



(pt) 
(GeV) 



) [pb/GeV] 



t (spin) 
,10600 — 
-7900 
,8600 



-2300 

,1620 

-1770 

,1870 

-900 

,1040 

-800 

,610 



118 
40 
52 
17.0 



Value 



bl < 0.75 
± (stat.) ± (syst.) 



± (spin) 



5.3 
5.8 

6.2 
6.7 
7.2 
7.7 
8.2 
8.7 
9.3 
9.7 
10.5 
11.5 
12.9 
14.9 
16.9 
19.7 
25.2 
33.9 
48.8 



26800 
19200 
13500 
12400 
8190 
6500 
4080 
3600 
2880 
2210 
1542 
1022 
531 
249 
119.2 
53.3 
15.9 
3.16 
0.407 



±5600 

±2800 

±1100 

±1100 

±610 

±400 

±280 

±200 

±140 

±100 

±51 

±35 

±16 

±10 

+6.7 

±3.0 

±1.1 

±0.43 

±0.084 



?:^ 



.0:0i2 

-0.017 



7.3 
7.8 
8.3 
8.8 
9.2 
9.8 
10.5 
11.5 
12.9 
14.9 
16.9 
19.7 
25.0 
33.6 
46.6 



9220 
7780 
4500 
3720 
3040 
2170 
1528 
1051 
528 
274 
136.2 
67.7 
16.9 
3.60 
0.462 



±980 

±720 

±320 

±270 

±280 

±140 

±59 

±39 

±17 

±12 

±7.5 

±3.6 

±1.4 

±0.48 

±0.093 



1140 

loSS 

|?8 



M 

-56 
+27 
-27 
+ 13.1 

l& 

Ms 

0.055 



5770 
2960 
3540 
2470 
1730 
1410 
1310 

m 

140 
'40 

fn 
;w 

-141 
,60 
-70 
.26.5 

:m 

% 
m 

0.049 
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Table 3: Non-prompt to inclusive production cross-section fraction fg as a function of Jjip pr for \y\j/^ < 0.75 under 
the assumption that prompt and non-prompt J/ip production is unpolarised (Ao = 0). The spin-alignment envelope spans 
the range of possible prompt cross-sections under various polarisation hypotheses, plus the range of non-prompt cross- 
sections within Ag = ±0.1. The first uncertainty is statistical, the second uncertainty is systematic, the third number is 
the uncertainty due to spin-ahgnment. 



Pt 
(GeV) 



6.0-7.0 
7.0-7.5 
7.5-8.0 
8.0-8.5 
8.5-9.0 
9.0-9.5 
9.5-10.0 
10.0-11.0 
11.0-12.0 
12.0-14.0 
14.0-16.0 
16.0-18.0 
18.0-22.0 
22.0-30.0 
30.0-70.0 



(Pt) 
(GeV) 



6.6 
7.3 
7.8 
8.3 
8.8 
9.2 
9.8 
10.5 
11.5 
12.9 
14.9 
16.9 
19.7 
25.0 
37.7 



Non-prompt to inclusive production fraction 
\y\ < 0.75 
/g + (stat.) + (syst.) ± (spin) 



0.175 
0.259 
0.236 
0.258 
0.291 
0.268 
0.320 
0.321 
0.327 
0.359 
0.405 
0.443 
0.479 
0.536 
0.656 



±0.057 
±0.038 
±0.030 
±0.032 
±0.030 
±0.025 
+0.026 
±0.018 
±0.019 
±0.017 
+0.024 
±0.030 
±0.030 
±0.039 
±0.059 



±0.032 
±0.002 
±0.007 
±0.017 
±0.005 
+0.008 
±0.006 
±0.007 
±0.003 
±0.003 
±0.008 
±0.005 
±0.004 
±0.008 
±0.008 



,0.064 
-0.062 
,0.066 
-0.080 
,0.061 
-0.076 
,0.054 
-0.074 
,0.058 
-0.079 
,0.054 
-0.076 
,0.062 
-0.083 
,0.050 
-0.077 
,0.051 
-0.078 
,0.044 
-0.069 
,0.046 
-0.072 
,0.048 
-0.073 
,0.040 
-0.063 
,0.032 
-0.050 
,0.030 
-0.045 



Table 4: Non-prompt to inclusive production cross-section fraction fs as a function of J/i/r pr for 0.75 < \y\j/^ < 1.5 
under the assumption that prompt and non-prompt J/)// production is unpolarised (-Id = 0). The spin-alignment envelope 
spans the range of possible prompt cross-sections under various polarisation hypotheses, plus the range of non-prompt 
cross-sections within Ag = ±0. 1. The first uncertainty is statistical, the second uncertainty is systematic, the third number 
is the uncertainty due to spin-alignment. 



Pt 



(Pt) 



Non-prompt to inclusive production fraction 
0.75 < |y| < 1.5 



(GeV) 


(GeV) 


fB 


± (stat.) 


± (syst.) 


± (spin) 


4.0-5.0 


4.7 


0.142 


±0.094 


±0.018 


,0.03!) 
-0.049 


5.0-5.5 


5.3 


0.183 


±0.049 


±0.036 


,0.039 
-0.058 


5.5-6.0 


5.8 


0.127 


±0.038 


±0.024 


,0.030 
-0.043 


6.0-6.5 


6.3 


0.188 


±0.033 


±0.019 


,0.042 
-0.057 


6.5-7.0 


6.8 


0.261 


±0.029 


±0.007 


,0.051 
-0.069 


7.0-7.5 


7.2 


0.230 


±0.025 


±0.017 


,0.041 
-0.061 


7.5-8.0 


7.8 


0.238 


±0.023 


±0.015 


,0.043 
-0.062 


8.0-8.5 


8.2 


0.226 


±0.022 


±0.032 


,0.036 
-0.055 


8.5-9.0 


8.8 


0.226 


+0.021 


±0.013 


,0.036 
-0.055 


9.0-9.5 


9.2 


0.261 


±0.021 


±0.009 


,0.040 
-0.060 


9.5-10.0 


9.8 


0.292 


±0.023 


±0.008 


,0.043 
-0.064 


10.0-11.0 


10.5 


0.315 


±0.016 


±0.004 


,0.040 
-0.061 


11.0-12.0 


11.5 


0.343 


±0.018 


±0.007 


,0.041 
-0.064 


12.0-14.0 


12.9 


0.352 


±0.016 


±0.005 


,0.033 
-0.054 


14.0-16.0 


14.9 


0.401 


±0.022 


±0.003 


,0.035 
-0.058 


16.0-18.0 


16.9 


0.450 


±0.031 


±0.006 


,0.036 
-0.058 


18.0-22.0 


19.7 


0.476 


±0.031 


±0.006 


,0.033 
-0.052 


22.0-30.0 


25.1 


0.542 


±0.042 


±0.015 


,0.029 
0.042 


30.0-70.0 


37.8 


0.594 


±0.060 


+0.016 


,0.029 
-0.040 
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Table 5: Non-prompt to inclusive production cross-section fraction fs as a function of J/t// px for 1.5 < lylj/,/, < 2 
under the assumption that prompt and non-prompt J/)/i production is unpolarised (Ae = 0). The spin-alignment envelope 
spans the range of possible prompt cross-sections under various polarisation hypotheses, plus the range of non-prompt 
cross-sections within = ±0.1. The first uncertainty is statistical, the second uncertainty is systematic, the third number 
is the uncertainty due to spin-alignment. 



Pt 



(Pt) 



Non-prompt to inclusive production fraction 
1.5 < |>i < 2 



(GeV) 


(GeV) 


fa 


+ (stat.) 


± (syst.) 


± (spin) 


1.0-4.0 


2.8 


0.100 


±0.053 


±0.039 


,0.061 
-0.031 


4.0-5.0 


4.6 


0.210 


±0.042 


±0.051 


,0.115 
-0.051 


5.0-5.5 


5.3 


0.218 


±0.043 


±0.006 


,0.097 
-0.050 


5.5-6.0 


5.8 


0.170 


±0.034 


±0.019 


,0.068 
-0.041 


6.0-6.5 


6.3 


0.180 


±0.034 


±0.048 


,0.057 
-0.042 


6.5-7.0 


6.8 


0.222 


±0.028 


±0.013 


,0.069 
-0.048 


7.0-7.5 


7.3 


0.195 


±0.025 


±0.017 


,0.049 
-0.044 


7.5-8.0 


7.8 


0.210 


±0.024 


±0.014 


,0.052 
-0.047 


8.0-8.5 


8.2 


0.216 


±0.022 


±0.022 


,0.042 
-0.044 


8.5-9.0 


8.8 


0.264 


±0.023 


±0.018 


,0.049 
-0.050 


9.0-9.5 


9.2 


0.287 


±0.026 


±0.015 


,0.051 
-0.052 


9.5-10.0 


9.7 


0.297 


±0.028 


±0.015 


,0.053 
-0.053 


10.0-11.0 


10.5 


0.335 


±0.019 


±0.004 


,0.043 
-0.055 


11.0-12.0 


11.5 


0.326 


±0.026 


±0.017 


,0.042 
-0.054 


12.0-14.0 


12.9 


0.357 


±0.022 


±0.015 


,0.034 
-0.045 


14.0-16.0 


14.9 


0.420 


±0.029 


±0.011 


,0.035 
-0.047 


16.0-18.0 


16.9 


0.517 


±0.038 


±0.007 


,0.039 
-0.048 


18.0-22.0 


19.7 


0.468 


±0.038 


±0.012 


,0.029 
-0.041 


22.0-30.0 


24.9 


0.605 


±0.058 


±0.005 


,0.021 
-0.032 



Table 6: Non-prompt to inclusive production cross-section fraction /g as a function of pr for 2 < lylj/,/, < 2.4 
under the assumption that prompt and non-prompt J/if> production is unpolarised (Ag = 0). The spin-alignment envelope 
spans the range of possible prompt cross-sections under various polarisation hypotheses, plus the range of non-prompt 
cross-sections within Ag = ±0.1. The first uncertainty is statistical, the second uncertainty is systematic, the third number 
is the imcertainty due to spin-aligimient. 



Pt 
(GeV) 



1.0-5.0 
5.0-6.0 
6.0-7.0 
7.0-7.5 
7.5-8.0 
8.0-8.5 
8.5-9.0 
9.0-9.5 
9.5-10.0 
10.0-11.0 
11.0-12.0 
12.0-14.0 
14.0-18.0 
18.0-30.0 



(Pt) 
(GeV) 



3.6 
5.5 

6.6 
7.2 
7.8 
8.2 
8.7 
9.2 
9.7 
10.5 
11.5 
12.9 
15.6 
21.7 



Non-prompt to inclusive production fraction 

2 < |>i < 2.4 
fs + (stat.) ± (syst.) ± (spin) 



0.098 
0.217 
0.289 
0.125 
0.231 
0.209 
0.183 
0.268 
0.249 
0.269 
0.297 
0.352 
0.348 
0.419 



±0.065 
±0.077 
±0.047 
±0.035 
±0.037 
±0.042 
±0.041 
±0.037 
±0.045 
±0.037 
±0.034 
±0.034 
±0.044 
±0.058 



±0.036 
±0.065 
±0.052 
±0.016 
±0.020 
±0.042 
±0.032 
±0.027 
±0.008 
±0.014 
±0.010 
±0.018 
±0.049 
±0.058 



,0.053 
-0.027 
,0.096 
-0.044 
,0.096 
-0.041 
,0.037 
-0.019 
,0.063 
-0.031 
,0.046 
-0.028 
,0.042 
-0.024 
,0.057 
-0.033 
,0.055 
-0.032 
,0.040 
-0.033 
,0.045 
-0.034 
,0.037 
-0.033 
,0.038 
-0.032 
,0.032 
-0.031 
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Table 7: Non-prompt Jlij/ production cross-sections as a function of J I if/ pr for \y\j/i/, < 0.75 under the assumption that 
prompt and non-prompt J/ifr production is unpolarised (/ig = 0), and the spin-alignment envelope spans the range of 
non-prompt cross-sections within = ±0.1. The first uncertainty is statistical, the second uncertainty is systematic. 
Comparison is made to FONLL predictions. 











^ p+ju-) [nb/GeV] 


Pt 


(Pt) 




\y\ < 0.75 






(GeV) 


(GeV) 


Value 


± (stat.) 


± (syst.) 


± (spin) 


FONLL prediction 


7.0-7.5 


7.3 


2.4 


±0.4 


±0.4 


±0.08 


1 7 + 0.7 
'-■ ' - 0.4 


7.5-8.0 


7.8 


1.8 


±0.3 


±0.3 


±0.06 


1 4 + 0.5 
^■^ - 0.3 


8.0-8.5 


8.3 


1.2 


±0.2 


±0.1 


±0.03 


1 1 + 0.4 
'^■^ - 0.3 


8.5-9.0 


8.8 


1.1 


±0.1 


±0.1 


±0.03 


9 + 0.3 

- 0.2 


9.0-9.5 


9.3 


0.8 


±0.1 


±0.1 


±0.02 


7 + "^ 

' - 0.2 


9.5-10.0 


9.8 


0.69 


±0.07 


±0.08 


±0.02 


62 + "-^2 

"■"^ - 0.15 


10.0-11.0 


10.5 


0.49 


±0.03 


±0.05 


±0.01 


47 + 

- 0.11 


11.0-12.0 


11.5 


0.34 


±0.02 


±0.04 


±0.01 


34 + 011 


12.0-14.0 


12.9 


0.19 


±0.01 


±0.02 


±0.004 


21 + 

"•^■^ - 0.05 


14.0-16.0 


14.9 


0.111 


±0.008 


±0.011 


±0.003 


1^ ^ ' - 0.024 


16.0-18.0 


16.9 


0.060 


±0.005 


±0.006 


±0.002 


069 + 


18.0-22.0 


19.7 


0.032 


+0.003 


±0.003 


±0.001 


035 + 


22.0-30.0 


25.0 


0.0091 


+0.0010 


+0.0011 


+0.0002 


0109 + " ""-^ 


30.0-70.0 


37.2 


0.0008 


±0.0001 


±0.0001 


±0.0000 


0009 + 



Table 8: Non-prompt J/ifr production cross-sections as a function of J/i/f pr for 0.75 < \y\jif < 1.5 under the assumption 
that prompt and non-prompt Jlij/ production is unpolarised (-Id = 0), and the spin-aUgnment envelope spans the range 
of non-prompt cross-sections within = ±0.1. The first uncertainty is statistical, the second uncertainty is systematic. 
Comparison is made to FONLL predictions. 











^ H*H-) [nb/GeV] 


Pt 


(Pt) 




0.75 < \y\ < 1.5 






(GeV) 


(GeV) 


Value 


± (stat.) 


± (syst.) 


± (spin) 


FONLL prediction 


5.0-5.5 


5.3 


4.9 


±1.7 


±1.2 


±0.15 


3.8 ±1? 


5.5-6.0 


5.8 


2.4 


±0.8 


±0.5 


±0.07 


3-0 ±i;^ 


6.0-6.5 


6.3 


2.5 


±0.5 


±0.4 


±0.07 


2.4 ± J:? 


6.5-7.0 


6.8 


3.3 


±0.5 


±0.5 


+0.09 


19 + 06 
^■^ - 0.5 


7.0-7.5 


7.2 


1.9 


±0.3 


±0.3 


±0.05 


1 5 + "-6 

- 04 


7.5-8.0 


7.8 


1.6 


±0.2 


±0.2 


±0.04 


'-■^ - 0.3 


8.0-8.5 


8.3 


0.9 


±0.1 


±0.1 


±0.02 


1-0 ±g;^ 


8.5-9.0 


8.8 


0.81 


±0.09 


±0.09 


±0.02 


0.83 ±«i« 


9.0-9.5 


9.3 


0.75 


±0.07 


±0.08 


±0.02 


0.68 ±0;?^ 


9.5-10.0 


9.8 


0.65 


±0.06 


±0.07 


±0.02 


56 + "-^o 


10.0-11.0 


10.5 


0.49 


±0.03 


±0.05 


±0.01 


43 + 


11.0-12.0 


11.5 


0.35 


±0.02 


±0.04 


±0.01 


30 + 


12.0-14.0 


12.9 


0.19 


±0.01 


±0.02 


±0.00 


19 + "o* 

"•^^ - 0.04 


14.0-16.0 


14.9 


0.100 


±0.007 


±0.011 


±0.002 


104 + 0029 


16.0-18.0 


16.9 


0.054 


±0.005 


±0.006 


±0.001 


061 +0016 


18.0-22.0 


19.7 


0.025 


±0.002 


±0.003 


±0.001 


030 + 

- 0.006 


22.0-30.0 


25.2 


0.0086 


±0.0009 


±0.0010 


±0.0001 


0093 + 


30.0-70.0 


38.0 


0.0007 


±0.0001 


±0.0001 


±0.0000 


0007 + 

' - 0.0001 
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Table 9: Non-prompt Jjiji production cross-sections as a function of Jjiff pj for 1.5 < < 2 under tlie assumption 
tliat prompt and non-prompt J/>// production is unpolarised (Ag = 0), and the spin-alignment envelope spans the range 
of non-prompt cross-sections within /ig = ±0.1. The first uncertainty is statistical, the second uncertainty is systematic. 
Comparison is made to FONLL predictions. 











^ //+/<-) [nb/GeV] 


Pt 


(Pt) 




1.5 < M < 2 






(GeV) 


(GeV) 


Value 


± (stat.) 


± (syst.) 


± (spin) 


FONLL prediction 


1.0-4.0 


2.8 


14.3 


±7.9 


±5.5 


+U.J 1 


Q /I _i_ 4.5 

9.4 ± 3 , 


4.0-5.0 


4.5 


8.3 


±2.0 


±1.9 


±0.21 


4 9 + 2.2 
^■^ - 1.6 


5.0-5.5 


5.2 


3.8 


±1.0 


±0.9 


±0.09 


3 4 + 1.5 

-'•^ - 1.0 


5.5-6.0 


5.7 


2.4 


±0.5 


±0.4 


±0.06 


^- ' - 0.7 


6.0-6.5 


6.3 


2.3 


±0.5 


±0.6 


±0.05 


2 1 + 0.9 
^- ^ - 0.6 


6.5-7.0 


6.8 


2.0 


±0.3 


±0.3 


±0.05 


1 7 + 

1 • ' - 0.4 


7.0-7.5 


7.3 


1.2 


±0.2 


±0.2 


±0.03 


^ 3 + 0.5 
^■^ ^ 0.4 


7.5-8.0 


7.8 


1.1 


±0.1 


±0.1 


±0.02 


1-1 


8.0-8.5 


8.3 


0.82 


±0.10 


±0.10 


±0.02 


0.87±«-33 


8.5-9.0 


8.8 


0.82 


±0.08 


±0.10 


±0.02 


0.71 ±«]^ 


9.0-9.5 


9.3 


0.65 


±0.07 


±0.07 


±0.01 


0.58 ±-i 


9.5-10.0 


9.8 


0.50 


±0.05 


±0.05 


±0.01 


0.48±«;|{ 


10.0-11.0 


10.5 


0.43 


±0.03 


±0.05 


±0.01 


36 + "'2 


11.0-12.0 


11.5 


0.25 


±0.02 


±0.03 


±0.01 


25 + 

"•^-^ - 0.06 


12.0-14.0 


12.9 


0.14 


±0.01 


±0.01 


±0.002 




14.0-16.0 
16.0-18.0 
18.0-22.0 
22.0-30.0 


14.9 
16.9 
19.6 
25.0 


0.081 
0.053 

0.023 
0.0064 


±0.007 
±0.005 

±0.002 
±0.0009 


±0.009 
±0.007 

±0.002 
±0.00081 


±0.001 
±0.001 

±0.0000 
±0.0001 


0.085 ±0:«?^ 
049 + 

024 + """'^ 
0071 +o'"»4 

U.WU/ 1 zr QQQJ2 



Table 10: Non-prompt Jjifi production cross-sections as a function of Jl^jj pr for 2 < lyly/^ < 2.4 under the assumption 
that prompt and non-prompt J/t// production is unpolarised (/le = 0), and the spin-alignment envelope spans the range 
of non-prompt cross-sections within Ae = ±0.1. The first uncertainty is statistical, the second uncertainty is systematic. 
Comparison is made to FONLL predictions. 



dpxdy 



Br(7/0- ^ [nb/GeV] 



Pt 


(Pt) 




2 < lyl < 2.4 




(GeV) 


(GeV) 


Value 


± (stat.) 


± (syst.) 


± (spin) 


5.0-6.0 


5.5 


3.2 


±1.3 


±1.0 


±0.04 


6.0-7.0 


6.5 


2.2 


±0.4 


±0.4 


±0.02 


7.0-7.5 


7.3 


0.5 


±0.2 


±0.1 


±0.01 


7.5-8.0 


7.8 


0.9 


±0.2 


±0.2 


±0.01 


8.0-8.5 


8.3 


0.6 


±0.1 


±0.1 


±0.01 


8.5-9.0 


8.8 


0.35 


±0.09 


±0.06 


±0.01 


9.0-9.5 


9.2 


0.39 


±0.06 


±0.06 


±0.01 


9.5-10.0 


9.7 


0.30 


±0.06 


±0.04 


±0.005 


10.0-11.0 


10.5 


0.22 


±0.03 


±0.03 


±0.003 


11.0-12.0 


11.5 


0.18 


±0.02 


±0.02 


±0.003 


12.0-14.0 


12.9 


0.11 


±0.01 


±0.01 


±0.002 


14.0-18.0 


15.6 


0.042 


±0.006 


±0.005 


±0.0004 


18.0-30.0 


21.3 


0.0059 


±0.0010 


±0.0007 


±0.0001 



FONLL prediction 



2.7: 



1.1 

0.7 

1 7 + 0.7 

1.2±?-5 



0.9: 



0.3 

0.4 
0.2 

7 + °' 

' - 0.2 

60 + "^^ 
49 + "i^ 

"•^^ - 0.12 

0.40 ± ? i* 



0.10 
0.10 
0.07 

21 + 

"•^-^ - 0.05 
0.04 
0.03 

053 + "-"'^ 
0.0097 ±!! S 



0.30: 



0.13: 
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Table 11: Prompt Jl^r production cross-sections as a function of Jlifi pj for lyU/,/, < 0.75. The central value assumes 
unpolarised (Ag = 0) prompt and non-prompt production, and the spin-alignment envelope spans the range of possi- 
ble prompt cross-sections under various polarisation hypotheses. The first quoted uncertainty is statistical, the second 
uncertainty is systematic. Comparison is made to the Colour Evaporation Model prediction. 









dpjdy 




ju+ju-) [nb/GeV] 


Pt 


(Pt) 




\y\ < 0.75 




(GeV) 


(GeV) 


Value 


± (stat.) 


± (syst.) 


± (spin) 


CEM prediction 


7.0-7.5 


7.3 


6.8 


±0.8 


+1.1 
-1.1 


+4.3 
-2.2 


2.8 


7.5-8.0 


7.8 


5.9 


±0.6 


+0.9 
-0.9 


+2.7 
-1.9 


2.2 


8.0-8.5 


8.3 


3.3 


±0.3 


+0.4 
-0.4 


+ 1.3 
-1.0 


1.7 


8.5-9.0 


8.8 


2.6 


±0.2 


+0.4 
-0.4 


+0.9 
-0.8 


1.3 


9.0-9.5 


9.2 


2.2 


±0.2 


+0.3 
-0.3 


+0.9 
-0.6 


1.0 


9.5-10.0 


9.8 


1.5 


±0.1 


+0.2 
-0.2 


+0.5 
-04 


0.8 


10.0-11.0 


10.5 


1.04 


±0.05 


+0.11 
-0.11 


+0.32 
-0.29 


0.60 


11.0-12.0 


11.5 


0.71 


±0.03 


+0.08 
-0.08 


+0.19 
-0.20 


0.41 


12.0-14.0 


12.9 


0.34 


±0.01 


+0.04 
-0.04 


+0.08 
-0.09 


0.24 


14.0-16.0 


14.9 


0.163 


±0.010 


+0.016 
-0.016 


+0.036 
-0.042 


0.128 


16.0-18.0 


16.9 


0.076 


±0.006 


+0.008 
-0.008 


+0.015 
-0.018 


0.071 


18.0-22.0 


19.7 


0.035 


+0.003 


+0.004 
-0.004 


+0.006 
-0.008 


0.035 


22.0-30.0 


25.0 


0.0078 


+0.0009 


+0.0009 
-0.0009 


+0.0010 
-0.0014 


0.0109 


30.0-70.0 


37.2 


0.0004 


±0.0001 


+0.0001 
-0.0001 


+0.0000 
-0.0000 


0.0008 



Table 12: Prompt Jjifi production cross-sections as a function of Jlifi pr for 0.75 < [yly/^ < 1.5. The central value 
assumes unpolarised (Zg = 0) prompt and non-prompt production, and the spin-alignment envelope spans the range 
of possible prompt cross-sections under various polarisation hypotheses. The first quoted uncertainty is statistical, the 
second uncertainty is systematic. Comparison is made to the Colour Evaporation Model prediction. 



Pt 
(GeV) 



5.0-5.5 
5.5-6.0 
6.0-6.5 
6.5-7.0 
7.0-7.5 
7.5-8.0 
8.0-8.5 
8.5-9.0 
9.0-9.5 
9.5-10.0 
10.0-11.0 
11.0-12.0 
12.0-14.0 
14.0-16.0 
16.0-18.0 
18.0-22.0 
22.0-30.0 
30.0-70.0 



(Pt) 
(GeV) 



5.3 
5.8 
6.2 
6.7 
7.2 
7.7 
8.2 
8.7 
9.3 
9.7 
10.5 
11.5 
12.9 
14.9 
16.9 
19.7 
25.2 
38.0 



dpxdy 



Br(7/0- ^ n*iJ.-) [nb/GeV] 
0.75 < \y\ < 1.5 



Value ± (stat.) ± (syst.) 



; (spin) 



21.9 

16.8 

11.0 

9.2 

6.3 

5.0 

3.2 

2.8 

2.1 

1.57 

1.06 

0.67 

0.34 

0.149 

0.066 

0.028 

0.0073 

0.0004 



±4.7 

±2.6 

±1.0 

±0.9 

±0.5 

±0.3 

±0.2 

±0.2 

±0.1 

±0.09 

±0.04 

±0.03 

±0.01 

±0.008 

±0.005 

±0.002 

±0.0008 

±0.0001 



+4.8 

-4.6 

+3.0 

-2.9 

+1.4 

-1.4 

+ 1.4 

-1.4 

+0.8 

-0.8 

+0.6 

-0.6 

+0.3 

-0.3 

+0.3 

-0.3 

+0.2 

-0.2 

+0.17 

-0.17 

+0.12 

-0.12 

+0.08 

-0.08 

+0.04 

-0.04 

+0.016 

-0.016 

+0.007 

-0.007 

+0.003 

-0.003 

+0.0008 

-0.0008 

+0.0001 

-0.0001 



^8.7 
-6.5 
^7.5 
-5.0 
^5.8 
-3.2 
^2.9 
-2.6 
^1.7 
-1.8 

(-1-2 
-14 

(-'■'* 

-0.7 

,0.8 

-0.6 

^0.5 

-0.5 

,0.30 

-0.35 

,0.19 

-0.24 

,0.12 

-0.15 

,0.06 

-0.08 

,0.024 

-0.031 

,0.009 

-0.013 

,0.004 

-0.005 

,0.0008 

-0.0011 

,0.0000 

-0.0000 



CEM prediction 



10.4 

7.2 

5.2 

3.7 

2.8 

2.1 

1.6 

1.3 

1.0 

0.79 

0.59 

0.39 

0.23 

0.120 

0.067 

0.032 

0.0100 

0.0007 
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Table 13: Prompt Jliji production cross-sections as a function of 7/(fr pr for 1.5 < |y|y/^ < 2. The central value 
assumes unpolarised (Ag = 0) prompt and non-prompt production, and tlie spin-alignment envelope spans tlie range 
of possible prompt cross-sections under various polarisation hypotheses. The first quoted uncertainty is statistical, the 
second uncertainty is systematic. Comparison is made to the Colour Evaporation Model prediction. 









dpfdy 


■Bi(J/i^ 


ju+ju-) [nb/GeV] 


Pt 


(Pt) 




1.5 < M < 2 




(GeV) 


(GeV) 


Value 


± (stat.) 


+ { svst ^ 


+ (spin) 


CEM prediction 


1.0-4.0 


2.8 


129 


±22 


+25 
-35 


+246 
-35 


43 


4.0-5.0 


4.5 


31.1 


±4.6 


+4.5 
-4.5 


+54.7 
-7.7 


17.7 


5.0-5.5 


5.2 


13.8 


±2.7 


_l_3.4 
-3.2 


+ 13.5 
-3.2 


10.0 


5.5-6.0 


5.7 


11.8 


±1.1 


±l:f 


+ 11.6 
-2.6 


6.7 


6.0-6.5 


6.3 


10.5 


±0.9 




+ 8.2 
-2.2 


4.8 


6.5-7.0 


6.8 


6.9 


±0.5 


+ 1.0 

-i.o 


+4.2 
-1.6 


3.4 


7.0-7.5 


7.2 


5.1 


±0.4 


+0.7 
-0.7 


+2.5 
-1.2 


2.6 


7.5-8.0 


7.7 


4.0 


±0.3 


+0.5 
-0.4 


+ 1.8 
-0.7 


1.9 


8.0-8.5 


8.3 


3.0 


±0.2 


+0.3 
-0.3 


+ 1.2 
-04 


1.5 


8.5-9.0 


8.7 


2.3 


±0.1 


+0.3 
-0.2 


+0.7 
-0.3 


1.2 


9.0-9.5 


9.2 


1.61 


±0.09 


+0.17 
-0.17 


+0.46 
-0.26 


0.89 


9.5-10.0 
10.0-11.0 


9.7 
10.5 


1.18 
0.86 


±0.08 
±0.04 


,0.13 
-0.12 

,0.10 
-0.09 


+0.32 
-0.21 
+0.21 
-0.16 


0.72 

0.53 


11.0-12.0 


11.5 


0.51 


±0.03 


+0.05 
-0.05 


+0.11 
-0.10 


0.35 


12.0-14.0 


12.9 


0.26 


±0.01 


+0.03 
-0.02 


+0.05 
-0.04 


0.21 


14.0-16.0 
16.0-18.0 

18.0-22.0 


14.9 
16.9 

19.6 


0.112 
0.050 

0.026 


±0.008 
±0.005 

±0.003 


+0.012 
-0.011 
+0.007 
-0.005 
+0.002 
-0.002 


+0.016 
-0.019 
+0.006 
-0.008 
+0.003 
-0.004 


0.106 
0.057 
0.028 


22.0-30.0 


25.0 


0.0042 


±0.0007 


+0.0005 
-0.0005 


+0.0003 
-0.0005 


0.0084 



Table 14: Prompt J/i// production cross-sections as a function of J/t/i pr for 2 < \y\j/if, < 2.4. The central value 
assumes unpolarised {/ig = 0) prompt and non-prompt production, and the spin-alignment envelope spans the range 
of possible prompt cross-sections under various polarisation hypotheses. The first quoted tmcertainty is statistical, the 
second uncertainty is systematic. Comparison is made to the Coloiff Evaporation Model prediction. 



Pt 



{Pt) 



^^•Br(7/^^/.V)[nb/GeV] 

2 < Ivl < 2.4 



(GeV) 


(GeV) 


Value 


± (stat.) 


± (syst.) 


± (spin) 


CEM prediction 


5.0-6.0 


5.5 


11.5 


±2.5 


+2.5 
-2.6 


+ 12.1 
-2.0 


7.8 


6.0-7.0 


6.5 


5.6 


±0.6 


+0.6 
-0.6 


+3.3 
-0.9 


3.9 


7.0-7.5 


7.3 


3.5 


±0.5 


+0.6 
-0.5 


+2.3 
-0.6 


2.3 


7.5-8.0 


7.7 


3.1 


±0.4 


+0.5 
-0.4 


+2.2 
-0.5 


1.8 


8.0-8.5 


8.3 


2.1 


±0.3 


+0.3 

-0.3 


+0.7 
-0.5 


1.4 


8.5-9.0 


8.7 


1.6 


±0.2 


+0.2 
-0.2 


+0.5 
-0.3 


1.1 


9.0-9.5 


9.2 


1.1 


±0.1 


+0.1 
-0.1 


+0.3 
-0.2 


0.9 


9.5-10.0 


9.7 


0.91 


±0.09 


+0.13 

-0.12 


+0.33 
-0.12 


0.68 


10.0-11.0 


10.5 


0.61 


±0.05 


+0.07 
-0.07 


+0.21 
-0.06 


0.47 


11.0-12.0 


11.5 


0.42 


±0.04 


+0.05 
-0.06 


+0.12 
-0.05 


0.32 


12.0-14.0 


12.9 


0.21 


±0.02 


,0.02 
-0.02 


+0.05 
-0.03 


0.18 


14.0-18.0 


15.6 


0.079 


±0.007 


+0.007 
-0.005 


+0.012 
-0.007 


0.071 


18.0-30.0 


21.2 


0.008 


±0.001 


+0.001 
-0.001 


+0.001 
-0.001 


0.012 
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